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Replication and meta-analysis of TMEM132D gene 
variants in panic disorder 

A Erhardt\ N Akula^, J Schumacher^, D Czamara\ N Karbalai\ B Muller-Myhsok\ 0 Mors"*, A Borglum^, AS Kristensen'*, 
DPD Woldbye^ P Koefoed^ E Eriksson^ E Maron^'2^ A Metspalu^'^°, J Nurnberger^\ RA Philibert^^ J Kennedy^^ K Domschke^^ 
A Reif^^ J Deckert^^ T Otowa^^ Y Kawamura^^ H Kaiya^^'^^ Y Okazaki^^ H Tanii^^, K Tokunaga2\ T Sasaki^^ JPA loannidis2^ 
FJ McMahon^ and EE Binder^ 

A recent genome-wide association study in patients with panic disorder (PD) identified a risk haplotype consisting of two single- 
nucleotide polymorphisms (SNPs) (rs7309727 and rs1 1060369) located in intron 3 of TMEM132D{o be associated with PD in three 
independent samples. Now we report a subsequent confirmation study using five additional PD case-control samples (n = 1670 
cases and n = 2266 controls) assembled as part of the Panic Disorder International Consortium (PanIC) study for a total of 2678 
cases and 3262 controls in the analysis. In the new independent samples of European ancestry (EA), the association of 
rs7309727 and the risk haplotype rs7309727-rs1 1060369 was, indeed, replicated, with the strongest signal coming from patients 
with primary PD, that is, patients without major psychiatric comorbidities (a7 = 1038 cases and /7 = 2411 controls). This finding 
was paralleled by the results of the meta-analysis across all samples, in which the risk haplotype and rs7309727 reached P-levels 
of P = 1 .4e - 8 and P = 1 .1 e - 8, respectively, when restricting the samples to individuals of EA with primary PD. In the Japanese 
sample no associations with PD could be found. The present results support the initial finding that TMEM132Dgene contributes 
to genetic susceptibility for PD in individuals of EA. Our results also indicate that patient ascertainment and genetic background 
could be important sources of heterogeneity modifying this association signal in different populations. 
Translational Psychiatry (2012) 2, e156; doi:10.1038/tp.2012.85; published online 4 September 2012 



Introduction 

TMEM132D was first identified as a possible candidate 
gene for panic disorder (PD) in a genome-wide association 
study conducted at the Max Planck Institute of Psychiatry 
in Munich. A two single-nucleotide polymorphism (SNP) 
haplotype in intron 3 of the gene containing the T allele 
of rs7309727 and the A allele of rsll 060369 was highly 
associated with PD in three German case-control samples, 
but did not reach genome-wide significance in the meta- 
analysis.'' Furthermore, three independent SNPs located 
in intron 4 were associated with the severity of anxiety 
symptoms in patients with a number of different primary 
psychiatric disorders. Additional evidence for involvement 
of the TMEM132D gene in anxiety comes from a family- 
based genome-wide study in patients with ADHD. Here, 
rsll 060369 was highly associated (P=2.3e-5) with 



symptoms of anxiety and depression in these patients as 
measured by the relevant subscales of the Child Behavior 
Checklist.^ Finally, recent data from our group also suggest 
that private non-synonymous coding mutations in TMEM132D 
are over-represented in patients with PD.^ These findings in 
humans supporting an involvement of TMEM132D in anxiety 
are corroborated by data in a mouse model of extremes in trait 
anxiety that showed a positive correlation of high anxiety- 
related behavior and expression of Tmem132d mRNA in 
the anterior cingulate cortex, as well as association of a 
Tmem132d SNP with anxiety-related behavior. The gene 
expression results in animals are consistent with the finding 
that the risk allele is correlated with higher TMEM132D 
expression in the frontal cortex in human postmortem brain 
tissue. These data would suggest that higher expression 
of TMEM132D is associated with increased measures of 
anxiety. 
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TMEM1 32D is a single-pass transinembrane protein, which 
is highly expressed in the cortical regions of the human and 
mouse brain. The function of the gene product is still unknown. 
It has been proposed that the protein may serve as a cell- 
surface marker for oligodendrocyte differentiation."^ Other 
recent data, however, show that it may be most prominently 
expressed in neurons and colocalizes with actin filaments,^ 
suggesting that TMEM132D may be implicated in neuronal 
sprouting and connectivity in brain regions important for 
anxiety-related behavior. 

As evidence from multiple studies and experimental 
approaches point to the putative relevance of TMEM132D in 
anxiety phenotypes, we wanted to test whether our initial 
association results could be replicated in independent cohorts 
for PD. We thus genotyped the top associated TMEM132D 
SNPs from Erhardt etal.^ in five additional PD samples drawn 
from the Panic Disorder International Consortium (PanlC).^ 
Finally, we performed a meta-analysis including all eight PD 
samples, three from the initial study and five additionally 
genotyped samples, for a total of 2678 PD cases and 3262 
controls. 



Materials and methods 

Samples. The four top associated case-control SNPs in 
the TMEM132D gene reported in the original study were 



genotyped in five additional PD samples (total number of 
PD patients n=1670 and controls n = 2266) as part of the 
Panic study.^ The contributing sites were: (1) Aarhus (cases 
n=^02, controls n = 511), Copenhagen (cases n=^4^, 
controls n = 345) and Gothenburg (cases n = 217) from 
Denmark and Sweden (ACG); (2) Tartu and Estonia (cases 
n = 217, controls n = 285), (3) Iowa (cases n = 38, controls 
n = 40) and Toronto (cases n=127, controls n=123) from 
USA/Canada (IT); (4) the National Institute of Mental Health 
(NIMH), USA (cases n = 77, controls n= 192); and (5) Japan 
(cases n = 760, controls n = 760) (Table 1). Subjects were 
matched for gender and age and must have had their first 
panic attack before developing an alcohol abuse disorder 
and before the age of 35. 

In the final meta-analysis, eight PD samples were included; 
three German samples are described in Erhardt et al? (MPI 
(discovery sample): recruited at the Max Planck Institute 
in Munich; PD German multicenter (Stage 2): recruited at 
Universities of Munster, Wurzburg, Bonn and Gottingen; Bonn 
(Stage 3): recruited at the University of Bonn). All patients in 
the eight PD samples had a DSM-IV or ICD-1 0 (Aarhus, MPP) 
based diagnosis of PD, with or without agoraphobia. Informa- 
tion regarding the recruitment and characterization of cases 
and controls is detailed in these references. ''"^"^ In most 
samples, PD was the primary diagnosis (ACG, Tartu, Toronto, 
Japan, MPP and PD German multicenter''), except for the 
cases from Iowa, the NIMH and Bonn. Cases and controls 



Table 1 Demographic data of eight PD samples included in the final meta-analysis 



Sample 


Number 


Gender 


PD diagnosis 


Comorbidity 




Cases 


Controls 


Cases 


Controls 


Cases 


Cases 


ACG total 
Aarhus 
Copenhagen 
Gothenburg 


460 
102 
141 
217 


856 
511 
345 
0 


Male 25.0% 
Female 75.0% 


Male 50.5% 
Female 49.5% 


PDA 56%, PD 44% 
PPD 77.6% 


MDE 22.2% 
SA 5.2% 
NS 47.1% 


Tartu 


217 


285 


Male 23.0% 
Female 77.0% 


Male 30.0% 
Female 70.0% 


PDA 61 .5%, PD 38.5% 
PPD 49.7 


MDE/BP 62.3% 
SA 15.6% 


IT total 
Iowa 
Toronto 


165 
38 
127 


163 
40 
123 


Male 29.7% 
Female 70.3% 


Male 26% 
Female 74% 


PDA 75.6%, PD 24.4% 
PPD 60.6% 


MDE/BP 39.9% 


NIMH 


77 


192 


Male 21.0% 
Female 79.0% 


Male 21.0% 
Female 79.0% 


PDA 38.8%, PD61.2% 
PPD 0% 


MDE/BP 100% 


Japan 


760 


770 


Male 31 .8% 
Female 68.2% 


Male 52.8% 
Female 47.2% 


PDA 57.3%, PD 42.7% 
PPD 35.8% 


MDE 42.8%, SA 5.2%, 
GAD 4.0%, GOD 3.3% 


MPI^ 


216 


222 


Male 33.2% 
Female 66.8% 


Male 32.9% 
Female 67.1% 


PDA 83.3%, PD 
16.7% PPD 100% 




PD German multicenter® 


225 


240 


Male 40.4% 
Female 59.6% 


Male 40.2% 
Female 59.8% 


PDA 67.5%, PD 
32.5% PPD 78.2% 


MDE 22.7% 

SA 8.9%, NS 25.3% 


Bonn® 


436 


452 


Male 33.3% 
Female 66.7% 


Male 50.9% 
Female 49.1% 


PDA 39%, PD 20.4% 
PPD 10.1% 


MDE/BP 73.4%, GAD 40.6% 
SchiD 10.6% 



Abbreviations: ACG, Aarhus, Copenhagen (Denmark), Gothenburg (Sweden); BP, bipolar disorder; IT, Iowa, Toronto (USA, Canada); MDE, depressive disorder; 
MPI, Max Planck Institute for Psychiatry, Munich; NIMH, National Institute for Mental Health (USA); NS, not specified; OAD, other anxiety disorders; PDA, panic 
disorder with agoraphobia; PD, panic disorder without agoraphobia; PD German multicenter, recruited in Munster, Wurzburg, Gottingen and Bonn; PPD, primary PD, 
exclusion of severe and recurrent major depression, bipolar disorder, schizophrenia and substance dependence; SA, substance abuse or dependence secondary to 
the panic disorder (mainly alcoholism); SchiD, schizophrenic disorder. 

Patients and controls were matched by age; mean age of all patients 39.1 (s.d. 1 1 .5); mean age controls 39.4 (s.d. 1 1 .9). 
^PD samples previously genotyped in Erhardt etal. 
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Table 2 Single-SNP association analysis, liaplotype analysis and meta-analysis results of the five newly genotyped samples from the PanIC 

Samples ACG Tartu IT NIMH Japan Combined 



Number of cases 
Number of controls 


460 
856 


217 
285 


165 
163 


77 
192 


760 
770 


1679 
2266 


SNPs 
rs7309727 

OR (CI) 


0.001 

1.36(1.13-1.65) 


0.5 

1 .09 (0.82-1 .47) 


0.7 

0.95 (0.66-1 .36) 


0.05 
0.63 (0.39-1.0) 


0.4 

0.94 (0.81-1.09) 


0.41 

1.04 (0.93-1.16) 


rs 11 060369 
P- value 
OR (CI) 


0.9 

0.99 (0.83-1.18) 


0.6 

0.92 (0.7-1.21) 


0.01 

0.67 (0.49-0.94) 


0.6 

0.90 (0.61-1.32) 


0.8 

0.98 (0.85-1.14) 


0.28 
0.94 (0.87-1.04) 


rs9804847 
P- value 
OR (CI) 


0.9 

1.0 (0.85-1.18) 


0.3 

1.15(0.88-1.49) 


0.05 
0.72 (0.52-1.01) 


0.17 
0.77 (0.53-1.11) 


0.2 

1.21 (0.93-1.58) 


0.97 
0.99 (0.89-1.11) 


rs4759958 
P- value 
OR (CI) 


0.07 
1.17(0.98-1.38) 


0.8 

0.96 (0.73-1 .28) 


0.1 

0.28 (0.93-1 .75) 


0.7 

1 .07 (0.73-1 .60 


0.3 

1 .07 (0.93-1 .24) 


0.036 

1.11 (1.0-1.23) 


Haplotype rs7309727-rs 11 060369 
Omnibus test 0.0054 


0.81 


0.21 


0.037 


0.33 


0.47 


CC 
P- value 
OR (CI) 


0.69 
0.98 (0.82-1.17) 


0.54 
0.92 (0.70-1.21) 


0.02 

0.67 (0.49-0.94) 


0.51 

0.99 (0.82-1.17) 


0.82 
0.98 (0.81-1.1) 


0.77 
0.99 (0.90-1.08) 


TA 
P- value 
OR (CI) 


0.0023 

1.36(1.12-1.65) 


0.5 

1.1 (0.81-1.49) 


0.7 

0.95 (0.66-0.37) 


0.057 
0.63 (0.4-1 .0) 


0.37 
0.93 (0.81-1.08) 


0.036 

1.11 (0.10-1.24) 


CA 
P- value 
OR (CI) 


0.006 

0.79 (0.67-0.93) 


0.9 

1.0 (0.99-1.01) 


0.01 

1.5(1.1-2.05) 


0.017 

1.6(1.1-2.39) 


0.16 
1.15(0.95-1.39) 


0.064 
0.91 (0.84-1 .0) 



Abbreviations: ACG, Aarhus, Copenhagen, Gothenburg; CI, 95% confidence interval; EA, European ancestry, all samples except the Japanese sample; IT, Iowa, 
Toronto (USA, Canada); NIMH, National Institute for Mental Health (USA); OR, odds ratio; PanIC, Panic Disorder International Consortium; SNP, single-nucleotide 
polymorphism. 

Nominally significant P-values are displayed in bold digits. 



from Iowa were recruited from a large longitudinal case and 
control adoption study of substance use. The samples from 
Bonn and NIMH were primarily recruited from families with 
bipolar disorder (Table 1). In the primary PD samples listed 
above, control individuals were recruited from the general 
population. The presence of PD and other psychiatric 
disorders was excluded by administering either a short 
questionnaire (Japan, ACG) or a structural interview (MPI: 
SCID, Toronto: Mini-SCID, Tartu: MINI). All studies were 
approved by the local ethical committees. 

Genotyping. The genotyping was performed at the NIMH. 
For genotyping, pre-designed TaqMan SNP assays were 
used according to the manufacturer's protocol on the MJ 
Research Tetrad PTC-225 Thermo Cycler (Bio-Rad Labora- 
tories, Hercules, CA, USA) and intensities were read on an 
Analyst HT fluorescent plate reader (Molecular Devices, 
Sunnyvale, OA, USA). The genotypes were scored using a 
clustering algorithm kindly provided by Sam Chen, Virginia 
Commonwealth University (Richmond, VA, USA). 

Statistics. Case-control association analyses were con- 
ducted using PLINK (http://pngu.mgh.harvard.edu/~purcell/ 
plink/) with an additive model for the single-SNP associations 
and the option hap-logistic for haplotype-based associations. 



P^O.05 were considered as nominally significant. Fixed- 
effects meta-analysis and forest plots were assessed 
with the R-library rmeta (http://www.r-project.org/). We first 
analyzed the newly genotyped five PD samples for single- 
SNP association and haplotype association (n= 1679 cases, 
n = 2266 controls). In the second step, a meta-analysis of 
these five samples was performed. Finally, a meta-analysis 
of all eight available PD samples, including the first three 
samples reported in Erhardt et al. (combined: cases 
n = 2678, controls n = 3262), was done, followed by meta- 
analyses restricted to European ancestry (EA) samples, 
only excluding the Japanese sample (EA: n=1863 cases, 
n = 24^^ controls). As the association of TMEM132D 
variants in the original publication was highest for PD without 
major comorbidities, we further stratified the EA samples by 
primary diagnosis of PD (PPD; cases n=1038, controls 
n = 2411), that is, excluding all patients with comorbid major 
depression, bipolar disorders, schizophrenia or substance 
dependence. 



Results 

Independent replication samples. Single-SNP analysis 
and haplotype analysis of the five newly genotyped samples 
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are presented in Table 2. Single-SNP associations in the 
same direction as in the initial three German samples were 
seen in the ACG sample and partly in the sample from IT for 
the top two SNPs, rs7309727 and rs11 060369, being 
strongest in the ACQ sample (for example, rs7309727, 
P= 0.001). The analysis of the haplotype rs7309727- 
rs11 060369 was nominally significant in the omnibus test 
in the ACG (P= 0.0054) and NIMH (P= 0.037) samples. The 
single-risk haplotype TA reported in the original study could 
be replicated with nominally significant results only in the 
ACG sample (P= 0.0023). In the meta-analysis of the five 
new PD samples, nominally significant results could be 
obtained only for rs4759958 (P= 0.036). Nominally signifi- 
cant associations in the same direction as in the first report 
were observed for rs7309727 for EA populations (P= 0.045) 
and PPD (P= 0.003), but not in the combined analysis 
(P=0.28). The meta-analysis of the haplotype rs7309727- 
rs11 060369 was not nominally significant using the omnibus 
test. However, the meta-analysis for the risk haplotype TA 
reached nominal significance (P= 0.036) with a trend for the 
protective haplotype OA (P=0.06), both with congruent 
directions to the first study. 

Meta-analysis across all samples. The results of the final 
meta-analysis including eight PD samples for single SNPs 
and haplotype rs7309727-rs 11 060369 are summarized in 
Table 3. 

In the meta-analysis of all EA samples with primary PD, the 
association with rs7309727 was strongest (P=1.8e-8, 
OR = 1.43) (Figure 1). Associations were less significant 
when analyzing all EA samples (P=2.5e-6) or all samples 
together, including the Japanese sample (P= 0.0007). 
When the meta-analysis was further restricted to primary 
PD samples recruited at European sites (cases n = 908, 
controls n = 2056), the significance increased to a P-value 
of 7.1e-9 (OR 1.48). Similar results were obtained for 
the risk haplotype TA for rs7309727-rs1 1060369 (EA PPD 
p=1.4e-8, OR 1.44; see Figure 2 and Table 3). It is 
noteworthy that the association of the protective haplotype 
CC is strengthened in all combined samples, especially in 
the EA sample, and that this haplotype is associated in the 
same direction with PD in all eight samples. 

Discussion 

The present results support the association of rs7309727 and 
the risk haplotype rs7309727-rs1 1060369 in the TMEM132D 
gene in EA samples selected for PPD, corroborating the 
initial finding that TMEM132D gene contributes to the genetic 
susceptibility for PD. 

In fact, the association of rs7309727 and the risk haplotype 
rs7309727-rs1 1060369 replicates in the new independent 
samples of European descent, with the strongest signal 
coming from patients with primary PD (see Table 2). This was 
paralleled by the findings from the meta-analysis across all 
samples in which the association reached a significance level 
below the threshold for multiple correction for genome-wide 
association studies in this subgroup of patients (see Table 3). 
This indicates that the association with the TMEM132D is 
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Figure 1 Combined meta-analysis for the SNP rs7309727 in the TMEM132D 
gene. ACG: Aarhus, Copenhagen, Gothenburg; Tartu (Estonia); European 
Ancestry: all samples except the Japanese sample; European Ancestry 
primary PD: exclusion of cases with severe and recurrent major depression, 
bipolar disorder, schizophrenia and substance dependence; IT: Iowa, Toronto 
(USA, Canada); MPI: Max Planck Institute for Psychiatry, Munich; NIMH: National 
Institute for Mental Health (USA). *Samples published in Erhardt etal.^ 
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Figure 2 Combined meta-analysis of the risk haplotype TA for rs7309727 and 
rsl 1060369. ACG: Aarhus, Copenhagen, Gothenburg; Tartu (Estonia); European 
Ancestry: all samples except the Japanese sample; European Ancestry primary PD: 
exclusion of cases with severe and recurrent major depression, bipolar disorder, 
schizophrenia and substance dependence; IT: Iowa, Toronto (USA, Canada); MPI: 
Max Planck Institute for Psychiatry, Munich; NIMH National Institute for Mental 
Health (USA). *Samples published in Erhardt etal.^ 



strongest when increasing the homogeneity of the samples 
with respect to both genetic background as well as case 
ascertainment and/or definition strategies. 

The association was not replicated in the studied 
Japanese sample, perhaps due to differences in risk alleles 
or population structure. Consistent with this idea, the tests 
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Table 3 Meta-analysis results for eight PD samples (PanIC: n = 5; Erhardt etal.: n = 3) for the top SNPs and the best associated haplotype rs7309727-rs1 1060369 in 
the TMEM132D gene 



SBmplos 


Combinod 


EuropoBn ancostry 


EuropoBn ancostry PPD 


Japan PPD 


Number of cases 


2678 


1863 


1038 


760 


Number of controls 


3262 


2411 


2411 


770 


SNPs 










rs7309727 










P- value 


0.0007 


3.98e — 6 


1 .05e — 8 


0.27 


OR (CI) 


1 16 M 05-1 26) 

I , i \J II ,\J\J 1 ,^\J J 


1 27 (^ 15-1 411 


1 45 M 20-1 721 


0 90 ^0 73-1 091 


p-het 


5 97e - 6 


0 0003 


NS 




rs 11 060369 










P- value 


0.0007 


7.64e -5 


0.0004 


0.78 


OR (CI) 


0.87 (0.74-0.94) 


0.83 (0.75-0.90) 


0.81 (0.73-0.89) 


0.97 (0.80-1.18) 


p-het 


0.01 


0.034 


NS 




rs9804847 










P- value 


0.05 


0.1 


0.12 


0.16 


OR (CI) 


1.08(1.0-1.18) 


1.07 (0.98-1.17) 


1.09 (0.96-1.22) 


1.29 (0.90-1.84) 


p-het 


0.009 


0.006 


0.005 




rs4759958 










P- value 


0.003 


0.006 


0.06 


0.08 


OR (CI) 


1.12 (1.1-1.21) 


1.14 (1.05-1.26) 


1.11 (0.97-1.26) 


1.20 (0.98-1.47) 


p-het 


NS 


NS 


NS 




Haplotype rs7309727-rs 11 060369 








Omnibus test 


0.00014 


1.12e-6 


5.02e - 8 


0.13 


CC 










P- value 


0.0003 


3.8e-5 


0.0004 


0.73 


OR (CI) 


0.87 (0.8-0.93) 


0.82 (0.75-0.90) 


0.80 (0.73-0.89) 


0.96 (0.80-1.17) 


p-het 


0.02 


0.022 


NS 




TA 










P- value 


0.001 


3.35e - 6 


1.4e-8 


0.24 


OR (CI) 


1.15(1.06-1.25) 


1.28(1.15-1.41) 


1.44(1.20-1.73) 


0.89 (0.74-1.06) 


p-het 


9.4e-6 


0.0002 


NS 




CA 










P- value 


0.9 


0.9 


0.09 


0.17 


OR (CI) 


1.0 (0.99-1.01) 


1.0 (0.99-1.01) 


0.91 (0.82-1.0) 


1.3 (0.92-1.83) 


p-het 


9.6e-4 


0.001 


NS 





Abbreviations: ACG. Aarhus, Copenhagen, Gothenburg; CI, 95% confidence interval; Combined, all eight PD samples including the initial German samples from 
Munich, PD German multicenter and Bonn; European ancestry, all samples except the Japanese sample; IT, Iowa, Toronto (USA, Canada); NIMH, National Institute 
for Mental Health (USA); NS, not significant; OR, odds ratio; p-het, P-value test for heterogeneity; PanIC, Panic Disorder International Consortium; PPD, primary 
panic disorder, exclusion of cases with severe and recurrent major depression, bipolar disorders, schizophrenia and substance dependence; SNP, single-nucleotide 
polymorphism. 

Results for the combined meta-analysis and subsamples. P-values surviving correction for multiple testing genome-wide are displayed in bold. 



of heterogeneity were significant in the combined sample, 
whereas they were not significant for the meta-analysis of the 
EA and European PPD samples. The differences in genetic 
background between the Asian and the European samples 
may contribute to this heterogeneity, as the functional variant 
underlying the association signal in intron 3 has not yet been 
identified. Further finemapping of the locus in both samples 
will be necessary to address this issue. Promising preliminary 
findings come from a GWAS for PD in a Japanese population, 
in which two SNPs in the TMEM132D gene located in the 
same intron as the top SNPs in the present study showed 
nominal association with PD."*^ This may suggest that a 
common but yet unidentified functional variant underlies the 
signals found in the EA and Japanese populations. It is 
increasingly recognized that significant differences in asso- 
ciations obtained from GWAS in different ancestry groups are 
notuncommon.^^ 

An interesting aspect of this study is that TMEM132D 
variants showed the strongest associations when restricting 



the European descent samples to patients with primary PD. 
While this is in line with our previous report, in which these 
variants were associated with PD but not with the dimensional 
measure of anxiety, it is at odds with a recent report that 
rs11 060369 is associated with dimensional measures of 
anxiety in patients with ADHD.^ In a recent pooled genome- 
wide association study in German individuals with personality 
disorders, one SNP in TMEM132D located in the intron 4 was 
nominally associated with Cluster B and C Personality 
Disorders (Weififlog et al., in revision). We had previously 
reported associations with dimensional measures of anxiety 
with a set of linked variants in intron 4 of TMEM132D\ha\ are 
not in linkage disequilibrium with the disease-associated 
variants in intron 3. In fact, these associations with dimen- 
sional anxiety measures were seen in patients with anxiety 
disorders, major depression as well as healthy individuals.'' 
Additional studies investigating different or common loci in 
TMEM132D associated with PD specifically or dimensional 
anxiety measures will be necessary. 
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This study presents genetic association data in the largest 
collection of patients with PD published to date (cases 
n = 2678, controls n=3262), within the PanIC consortium. 
While we observe interesting supportive data for the asso- 
ciation of TMEM132D with primary PD in individuals of 
European ancestry, additional studies in larger samples with 
an increased marker density are needed to identify the 
possible sources of heterogeneity of the observed associa- 
tion, and seek further replication and extensions of our 
findings in samples of different ethnic backgrounds, and to 
identify the putative functional variants. 
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